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ABSTRACT 



We report on a 63-ks long XMM-Newton observation of the accreting millisecond pulsar SAX J1808. 4-3658 during the 
latest X-ray outburst which started on September 21st 2008. The pn spectrum shows a highly significant emission line 
in the energy band where the iron K-a line is expected, and which we identify as emission from neutral (or mildly 
ionized) iron. The line profile appears to be quite broad (more than 1 keV FWHM) and asymmetric; the most probable 
explanation for this profile is Doppler and relativistic broadening from the inner accretion disc. From a fit with a diskline 
profile we find an inner radius of the disc of 8. 7^7 R-s> corresponding to 18.0^g'g km for a 1.4 M neutron star. The 
disc therefore appears truncated inside the corotation radius (31 km for SAX J1808. 4-3658) in agreement with the fact 
that the source was still showing pulsations during the XMM-Newton observation. 
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1. Introduction 

Accreting millisecond pulsars (AMSP) are generally in- 
terpreted as the evolutionary link between the low mass 
X-ray binaries (LMXB) and rotation powered millisec- 
ond pulsars. According to the recycling scenario, the 
latter are in fact formed after a phase of mass accre- 
tion from a low mass companion star, which ultimately 
spins the neutron star (NS) up to m s periods (see e.g. 
iBhattacharva fc van den Heuve]|fT99ll). SAX J1808.4-3658 



was t he first AMSP discovered (jWiinands fc van der Klid 
1998). In ten years the population of AMSP has grown 
to ten sources, all residing in close and transient binaries, 
■ but SAX J1808. 4-3658 can still be considered the corner- 
stone of its class as it has repeatedly gone into outburst 
almost every two years since its first detection, making it 
the most observationally rich source of its class. Its tim- 
ing and orbital properties have bee n extensively studie d 
and are still de b ated t o some extent (iBurderi et a l. ( 2006), 



iDi Salvo et "all (|2008D . lHartman et alj (|2008D .~H08 here 



inafter). Its broadband X-ray spectral emission during the 
outburst phase wa s studied using RXTE observations of 
the 199 8 outburst (iGilfanov et al.l (|1998[ ) , iHeindl fc Srmthl 
(|1998|) . lGierlinski et alj 1)20020 . GDB hereinafter). 

GDB decomposed the 3.0-150.0 keV X-ray spectrum 
into a soft component, originating either from the disc or 
the NS surface, and a dominating hard component extend- 
ing to energies up to 100 keV (kT e > 35 keV), interpreted as 
Comptonization of soft photons, probably in the accretion 
column above the magnetic caps of the NS. A disc reflection 



component was also observed producing the typical hump 
at 30 keV and a Ka iron line at ~ 6.5 keV. 

Broad iron emission lines, in the 6.4-6.97 keV en- 
ergy range, ha ve been obs erved from many NS LMXB 
(iDi Salvo et alJ (120051) . iBhattacharvva fc Strohmave: 



(I2007D. ICackett et al.l (|2008bD . iPandel et all ([200 



D'A i et alj ( 20081 ). and references therein). Their broaden- 
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ing and in particular the asymmetry of the line is generally 
interpreted as the effect of the relativistic Keplerian 
motion of the reflecting plasma in an optically thick, and 
geometrically thin, accretion disc, i n the deep gravitat ional 
well of the nearby compact object (jFabian et al.lll989l) . 

The determination of the physical properties that pro- 
duce the broadening of the iron line in an X-ray pulsar 
like SAX J1808. 4-3658 is invaluable for accretion theories 
of fast rotators, as it gives constraints on one of the key 
parameters, the radius at which the disc terminates as the 
magnetosphere starts to lift off matter towards the mag- 
netic poles. Motivated by this, we obtained a 63 ks XMM- 
Newton Target of Opportunity (ToO) observation of this 
source, during the latest outburst in September 2008. 



2. Observation and data reduction 

SAX J1808.4-3658 was found in outburst on 2008 
September 21 by RXTE , and since then it has been ob- 
ject of an intensive observational campaign. A preliminary 
analysis of the 2-10 keV Swift XRT publicly available light 
curve shows that the outburst had its maximum around 
September 24. 

XMM-Newton observed SAX J1808.4-3658 as a ToO 
observation for 63 ks on 2008 October 1 (start time MJD 
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54739.99517), roughly one week after the assumed outburst 
peak. The EPIC-pn camera operated in timing mode, to 
prevent photon pile-up and to allow an analysis of the co- 
herent and aperiodic timing behaviour of the source. The 
same observing mode was used for EPIC/MOS2 CCDs, 
while EPIC/MOS1 was operated in small window mode, 
and the RGS in spectroscopy mode. A bright external flare 
in the background is present during the first 2 ks of the ob- 
servation, so that we excluded this time interval from our 
analysis. 

The high source flux saturated the XMM-Newton 
telemetry rate, so that the XMM-Newton Science 
Operations Centre decided to switch off one of the MOS 
cameras (MOS2) in order to allocate more band to the 
EPIC-pn instrument, roughly 35 ks after the beginning of 
the observation. The EPIC-pn spectrum showed no notice- 
able difference before and after the MOS-2 turning off, and 
the whole EPIC-pn data set is therefore considered for our 
analysis. The MOS spectra were severely affected by pile 
up, and we did not consider them for our analysis. 

Data were extracted and reduced using SAS v. 8. 0.0. 
We produced a calibrated EPIC-pn event list through 
the epproc pipeline, using the most updated calibration 
tools (epfast tool). We extracted the source spectrum se- 
lecting a 13 pixel wide stripe around the source position 
(RAWX=37) equivalent to 53" .3 (which should encircle 
more than 90% of the energy up to 9 keVQ), and consid- 
ered only PATTERN < 4 and FLAG=0 events. The back- 
ground spectrum was extracted in the RAWX=4-15 and 
RAWY=2-198 CCD region. 

The EPIC-pn energy channels were grouped with a 
compression factor of three so as not to oversample the in- 
strument energy resolution. The count rate observed in the 
EPIC-pn CCD showed a slight increase during the obser- 
vation from ~ 700 c/s, at the beginning of the observation, 
to 780 c/s at the end of the observation. 

RGS data were processed using the rsgproc pipeline to 
produce calibrated event lists, first and second order spectra 
and response matrices. The net count rates observed by 
the RGS1 and RGS2 in their first order, which are the only 
considered here, are 17.8 c/s and 20.0 c/s respectively RGS 
spectra were binned in order to have at least 25 counts/bin. 

To have a benchmark of the pn response, we also 
analysed, using the standard pipelines, a 1.1 ks Swift 
XRT observation taken in window timing mode (Obsld 
00030034033), which started on MJD 54740.60118, and is 
therefore simultaneous with the XMM-Newton pointing. 
The X-ray spectral package we use to model the observed 
emission is HEASARC XSPEC v. 12.4. 

3. Spectral analysis and results 

In our spectral analysis, we have first considered the com- 
bined spectrum of the RGS 1-2 (0.6-2.0 keV) and of the 
EPIC-pn (2.0-11.0 keV) data. The choice of the lower 
bound of the RGS is due to the presence of several cal- 
ibration residuals at ~ 0.5 keV we chose not to model. 
The model we have first focused on is composed of a multi 
temperature disc emission (diskbb) and a power law. We 
model the effect of the interstellar photoelectric absorp- 
tion using the phabs component. The addition of another 

1 see XMM-Newton Users handbook, issue 2.6, available at 
|http://xmm.esac.esa.int/external/xmm_user_support / 



Table 1. Fitting parameters of the 0.6-11.0 keV combined 
RGS + EPIC-pn spectrum (centre) and of the 1.4-11.0 
keV EPIC-pn spectrum (right) of SAX J1808.4-3658. The 
continuum was modelled by a diskbb (DBB), a bbodyrad 
(BB) and a power law (PL). Errors on each parameter are 
quoted at the 90% confidence level, as for the rest of the 
paper. 
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" The absorption column in the EPIC-pn spectrum was varied 
in the range indicated by the RGS-EPIC-pn spectrum. 
6 Absorbed flux in the 2.0-10.0 keV band 
c Unabsorbed flux extrapolated to the 0.05-150 keV band 



(single temperature) soft component is strongly required 
by the data, as the chi squared improves by A% 2 = 1104 
with the addition of just two parameters for 3156 degrees 
of freedom (d.o.f.). An edge at 0.871 keV, identified with 
O VIII absorption, is present in the RGS data at a high 
significance (Ax 2 = 160 for the addition of 1 parameter). 
A line at 6.5 keV also clearly emerges in the residuals (see 
FigO}, at the energy of Ka fluorescent emission of neutral 
or mildly ionized iron. We obtain a A\ 2 = 399 for the ad- 
dition of 3 parameters when we include a Gaussian to fit 
this component. The line is indeed very broad, as we ob- 
tain a = 1.1 ± 0.2 keV. This broadness, together with the 
apparent asymmetric profile, leads us to consider a model 
in which the iron emission is broadened by relativistic mo- 
tion near the NS. W e therefore find that a diskline model 
(Fab ian et al]|1989f ) fits this feature well, giving an inner 
disc radius consistent with the range expected from accre- 
tion theories onto a pulsar (6 — 15 R g in the case of SAX 
J1808. 4-3658, see discussion). This model, together with 
residuals, is plotted in FigJ^Jand the obtained best fit pa- 
rameters are listed in the middle column of Table [TJ 

The final chi squared of the fit in the 0.6-11.0 keV is 
still quite large and this is due to the presence of several 
features, especially in the RGS data, which are probably 
of calibration origin (e.g. at 0.7 and 0.85 keV, see residu- 
als of Fig(2]) , and which we choose not to fit. Considering 
that the main goal of the observation (and of this Letter) is 
the study of the iron Ka line, in order to check the stabil- 
ity of the obtained model, and to prevent spurious effects 
from a possible wrong cross-calibration between the instru- 
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Fig. 1. Residuals, in units of a, of the 1.4-11.0 keV EPIC- 
pn spectrum, with respect to the best fit continuum model 
(right column of Table [T]) when the iron line is not included 
in the fit. 

ments, we have also analysed the EPIC-pn dataset alone, 
in a bandwidth extending to the lowest energy allowed. We 
are forced to limit the energy range to the 1.4-11.0 keV 
band, as below large residuals (> 6tr) appear with respect 
to any reasonable continuum model, as we could check with 
a simultaneous fit with the Swift XRT dataset. We fit two 
evident absorption features at 1.8 and 2.2 keV, which are 
probably of instrumental origin as they lie near the Si K and 
Au M edges, which are known to influence the effective area 
determination if incorrectly calibrated. However, we take 
advantage of the low energy spectral information given by 
the RGS data, by taking the absorption column to vary in 
the range determined by the best fit in the whole RGS/pn 
energy range. A quick glance at Table [T] confirms how the 
best fit parameters do not change significantly fitting the pn 
alone. This is especially true for the hard power-law com- 
ponent and the iron emission line parameters, which arc 
indeed perfectly consistent with those obtained from the fit 
in the 0.6-11 keV band. 

We note that the diskline model gives a slightly bet- 
ter fit for the Fe-Ka line than a symmetric Gaussian pro- 
file, with an improvement of A\ 2 = 13 for the addition 
of three parameters (over 588 d.o.f.) with respect to the 
Gaussian modelling. Together with the evident difficulties 
in explaining the observed broadening of this feature with 
alternative models (see discussion), this definitely leads us 
to interpret the line profile as being produced in the inner 
regions of the accretion disc. To be conservative we discuss 
our results in the light of both the best fits presented in 
Table 1, which we remark differ significantly only for the 
soft components, which themselves are less constrained by 
the EPIC-pn alone. 

To obtain an estimate of the unabsorbed bolometric flux 
we extrapolate the model of Table [1] to the 0.05-150 keV 
band. We are confident about our estimate, as values of 
a ~ 2 are easily observed in SAX J1808. 4-3658 (Gilfanov et 
al. 1998). We finally remark that the 2-10 keV flux estimate 
we obtain is in good accord with the one that can be derived 
by fitting our model to the Swift XRT spectrum (Fxrt(2 — 
10) = (7.8 ± 0.2) x 10~ 10 erg cm" 2 s" 1 ). 
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Fig. 2. 0.6-11.0 keV XMM-Newton combined RGS+EPIC- 
pn spectrum of SAX J1808. 4-3658 (top) and residuals, in 
units of er, with respect to the best fit model given in the 
central column of Table [T] Blue refers to the EPIC-pn, red 
to RGS1 and green to RGS2 data and residuals. The model 
components are also plotted: the power law as a magenta 
solid line, diskbb as a light blue dashed line, bbodyrad as 
an orange dash-dotted line and the Ka diskline compo- 
nent as a dotted magenta line. 

Once the orbital motion is corrected (|Patruno et alj 

120081 ). the well known coherent pulsation is easily detected 
in the pn data at a frequency v = 400.9752102 Hz, con- 
sistent with the 2005 value (H08). The power spectrum at 
frequencies lower than 300 Hz is presented in Fig|31 and is in 
line wit h the typical power spectr a observed for this source 
(see e.g. Ivan Straaten et al.ll2003f ). Three Lorentzians, with 
central frequencies of 0.92(5) Hz, 1.56(18) Hz, and 34.7(6) 
Hz, respectively, are used to model the timing features we 
observe. The first two (labelled A and B in the figure) de- 
scribe the broad bump already observed below 10 Hz dur- 
ing other outbursts of this source. The third Lorentzian, la- 
belled C, could be tentatively identified with the component 
observ ed during the 2002 outburst by Ivan Straaten et all 
(2005]), who named it L^. 

4. Discussion 

The spectrum of the AMSP SAX J1808.4-3658 observed 
by XMM-Newton is of extremely good quality. For the first 
time this opens the possibility of exploiting its high sensitiv- 
ity and spectral resolution capabilities to study the inner- 
most regions of the accretion disc around a rapidly rotating 
pulsar. The estimate of the inner disc radius we obtain from 
the broadening of the iron Ka fluorescent line is invaluable 
to test the theories of accretion onto a fast object. This can 
be used to constrain the truncation disc radius, which is ex- 
pected to lie within the centrifugal boundary represented by 
the corotation radius, as indicated by the presence of pulsa- 
tions at the known 401 Hz frequency in the XMM-Newton 
light curve. Moreover it will be useful for an estimate of the 
lever arm of the torque exerted on the NS because of accre- 
tion, when a timing solution covering the whole outburst 
will be available. 

Despite that the measured asymmetry of the iron line 
does not exclude a Gaussian profile on a statistical basis, we 
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Fig. 3. (Top panel) White-noise subtracted PSD of the 63 
ks XMM-Newton observation of SAX J1808.4-3658. The 
time series was sampled at 125 fis corresponding to a 
Nyquist frequency of 4000 Hz. The PSD was obtained aver- 
aging about 500 power spectra of 2 x 10 20 temporal bins (~ 
133 s) each. The white-noise (1.97423(25)) was estimated 
by fitting the spectrum between 2 and 3 kHz with a con- 
stant. The spectrum was logarithmically rebinned using a 
factor 1.01. Spurious frequencies related to the 1 67 Hz in- 
strumental signal (see lJuett fe Chakrabartvll2003T ) were re- 
moved from the spectrum. The three components labelled 
as A, B and C are the Lorentzians used to fit the spectrum. 
(Bottom panel) Residuals in units of a. 



propose an interpretation of its large broadening in terms 
of the relativistic motion of the reflecting material in the 
space-time bent by the gravitational influence of the NS. 
Disc reflection in fact seems the only viable option to ex- 
plain the large broadening and the low ionization state mea- 
sured for the emitting material, in a source whose broad- 
band spectrum is dominated by a hard (high temperature, 
low optical depth) Comptonized component (GDB). 

For X-ray pulsars, the inner disc radius, Rd, has to 
meet some observational constraints implied by the pres- 
ence of pulsations. Rd must be larger than the NS sur- 
face, i?„ ~ 10 km depending on the equation of state, and 
not much larger than the corotation radius, defined as the 
point where the magnetosphere rotation equals that of an 
assumed Keplerian disc, R c = (G MP 2 /Air 2 ) 1 / 3 = 3lm\^ 
km for SAX J1808. 4-3658, where m L4 is the NS mass in 
units of 1.4 M©. If not met, the latter condition implies 
that the magnetic field is not able to capture matter be- 
cause it spins more rapidly than the disc, therefore result- 
ing in a centrifugal bar rier. Theoretical underst anding and 
3D MHD simulations (R omanova et al.l (|2008f ) and refer- 
ences therein) showed how the corotation condition is not 
to be taken strictly, as the magnetic field is able to re- 
lease enough angular momentum to the disc to make it 
leap over the centrifugal barrier. In any case an inner disc 
radius well beyond the corotation radius would imply an 
efficient propeller ejection of matter, which seems not to be 



the case for SAX J1808. 4-3658 (which is observed to show 
pulsations during this observation, and therefore to accrete 
matter onto the magnetic poles). 

The estimate we made for the inner disc radius from the 
iron line broadening is Ri n = 18.0^g'g mi. 4 km. While the 
lower bound essentially overlaps the NS radius, the upper 
limit we give at a 90% confidence level as Rt n < 25.6 mi. 4 
km is in agreement with the expectation that the inner disc 
radius is inside the corotation boundary. Also in the case 
of a more massive 2.0 M Q NS, our upper limit is still of the 

order of R c , R^ ^ /R^' ^ < 1.05. This measurement thus 
well fits in the small zone available for a fast rotating NS 
to efficiently accrete matter. 

The line modelling holds if we consider both the broader 
0.6-11.0 keV RGS+EPIC-pn spectrum and the restricted 
1.4-11 keV bandwidth. The spectral information in the for- 
mer is nevertheless fundamental to constrain the absorption 
column to SAX J1808. 4-3658. The value we obtain is signif- 
icantly larger than the average Galactic value (0.13 x 10 22 
cm -2 ). This could in principle indicate the presence of ad- 
ditional neutral absorbers in the proximity of the source, 
possibly revealing ma ss lost from the syst em during the 
Roche Lobe over flow dDi Salvo et al.ll2008t) . The presence 
of an evident O VIII edge in the RGS data also indicates 
the presence of ionized absorbing matter along the line of 
sight. 

The measured continuum is in good agreement with 
the expectations for these kinds of sources (GDB, 
iGierlihski fc Poutanenl (|2005h ). We clearly detect two soft 
components, interpreted as coming from the accretion disc 
and from the NS surface (or a good fraction of it), as well as 
a dominant hard component, which is well described by a 
simple power-law in the XMM-Newton energy band. To de- 
rive the bolometric X-ray luminosity, we use our estimate 
of the 0.05-150.0 keV bolometric flux (see Table [T]), and 
obtain Lx = 6.6 x 10 36 d 2 5 , where d3.5 is the distance t o 
the source in units of 3.5 kpc (|Gallowav fc Cummind l2006') . 
This estimate is in line with the bolometri c flux measured 
for pre vious outbursts of this source (see e.g. lHartman et ail 
()2008D ). In the hypothesis the X-ray luminosity well tracks 
the mass accretion rate onto the NS, we obtain a mass ac- 
cretion rate of M ~ 5.6 x 10~ 10 d§ 5 R 6 m^ M Q /yr at 
this stage of the outburst, an estimate however affected by 
the uncertainties implied by its derivation from a tighter 
bandwidth spectrum. 

The parameters measured for the soft components (see 
Table [T]) can be used to obtain further constraints on the 
geometry of the emitting regions. The normalization of 
the multi temperature accretion disc black body (diskbb), 
Ndbb, can be used to obtain an estimate of the apparent 
inner radius r;„ = 0.35\/Ndbb I cos7 d.3,5 km, which is in- 
deed of the order of Rt n (see e.g. iKubota et al.l (|1998( )). In 
our case we obtain different values from the fit in the two 
considered spectral bands, ri n ~40-70 (cosi) -1 / 2 km, and 
large uncertainties are still present in the value of the incli- 
nation angle i of the system (i J> 60° as derived from the 
disklinc fit). Despite still being affected by large uncertain- 
ties, the range of values we derive for the inner disc radius 
from the diskbb component is still quite reasonable, and in 
agreement with the measure derived from the line fitting. 

The single temperature blackbody component we detect 
with a high statistical significance in our spectral fits most 
probably comes from the NS surface (or a good fraction 
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of it) as indicated by the radius of the blackbody emitting Wijnands, R. & van der Klis, M. 1998, Nature, 394, 344 

region which is 6 — 10 km (for a distance to the source of 3.5 
kpc, see Table 1), in line with the expectations of a hotspot 
origin. 

Thus, in this Letter we have reported the results 
obtained from very high quality XMM-Newton spectra. 
Although some systematic features are still present in these 
spectra (due to the high statistics), the best fit model is sta- 
ble and the best fit parameters do not dramatically depend 
on the energy band (or instrument) in which the spectral 
fit is performed. The most important result we report here 
concerns the observation of a broad emission line at ~ 6.5 
keV. The most probable interpretation of this feature is 
reflection of the hard Comptonization spectrum in the in- 
nermost accretion disc. The inner disc radius derived in 
this way confirms the expectation that the accretion disc 
is truncated inside the corotation radius. This is the first 
time such an estimate can be made for a fast rotating X- 
ray pulsar, and confirms the basic elements of theories of 
accretion onto these objects. 

During the preparation of this manuscript we became 
aware that a ~ 40 ks Suzaku obs ervation also detecte d 
the iron line object of this paper (jCackett et alj [200 8a) . 
Although the parameters they found for the continuum are 
slightly different to those reported in this Letter, the con- 
clusions stated there about the line modelling confirm our 
findings. 
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